Passive swimming in low Reynolds number flows 
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The possibility of microscopic swimming by extraction of energy from an external flow is discussed, 
focusing on the migration of a simple trimer across a linear shear flow. The geometric properties 
of swimming, together with the possible generalization to the case of a vesicle, are analyzed. The 
mechanism of energy extraction from the flow appears to be the generalization to a discrete swimmer 
of the tank-treading regime of a vesicle. The swimmer takes advantage of the external flow by both 
extracting energy for swimming and "sailing" through it. The migration velocity is found to scale 
linearly in the stroke amplitude, and not quadratically as in a quiescent fluid. This effect turns out 
to be connected with the non-applicability of the scallop theorem in the presence of external flow 
fields. 
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One of the problems that will have to be solved in 
the realization of microscopic artificial swimmers ("mi- 
crobots") is the energy source for locomotion. One pos- 
sibility, is that of an external drive, e.g. by means of 
micromagnets inserted in the device [1]. The final goal, 
however, is that of autonomous microbots that are able 
to migrate, like their biologic companions, in response to 
external stimuli such as gradients in a chemical or tem- 
perature field. 

Many situations involve swimming in a non-quiescent 
fluid, say the interior of a blood vessel or the neighbor- 
hood of other swimming microorganisms. These situa- 
tions involve the presence of gradients in the fluid ve- 
locity and the possibility, in principle, of energy conver- 
sion from external flow to swimming. A similar idea has 
been put forward as regards the possibility of microscopic 
swimming aided by rectification of thermal fluctuations 
[2]. 

Of course, unless imagining some sort of energy storage 
in the microbot, the swimming velocity would be quite 
small, of the order of the fluid velocity difference across 
the device body. Nevertheless, this may be not an issue 
if long time or collective behaviors are sought. Examples 
in nature of such behaviors do exist. A most notable one 
is the Fahraeus-Lindwist effect: red blood cells in small 
arteries passively adapt their shape to the flow and align 
in the middle of the vessel, thus decreasing the vessel 
fluid- mechanic resistivity [3]. 

Microscopic swimming takes places in an environment 
dominated by viscous stresses, in which geometrical as- 
pects play a dominant role. A strategy that has been 
utilized to highlight these aspects, has been to focus on 
systems with a minimal number of degrees of freedom. 
Beyond various kinds of squirmers [4], examples are the 
three-sphere swimmer [5, 6], the " pushmepullyou" [7] and 
PurccU's three-link swimmer [8, 9]. 

It is an easy guess that some of the geometrical con- 
strains valid in a quiescent fluid, will cease to hold in 
the presence of an external flow. In particular, the scal- 
lop theorem statement [8], that reversing a sequence of 
deformations brings back the swimmer to its original po- 



FIG. 1: The coordinate system. The bead coordi- 
nates of the undeformed trimer in the rotating frame, are 
^[°'/R = (1/V3, 0, 0), x™/i? = -(1/(2^3), 1/2, 0), //? = 
(-1/(2V^), 1/2,0). 



sition, will cease to be true. 

We are going to discuss an example of device that 
swims extracting energy out of an external flow, based 
on a generalization of the three-spheres swimmer of [5] 
(see Fig. 1). The swimmer is placed in an unbounded 
linear shear flow 



U(X) = U(0) + (0,aXi,0), 



(1) 



and wants to migrate along the gradient direction Xi. 
We assume absence of external forces and torques on the 
trimer, except for those by the fluid, and neglect the 
interaction between the trimer arms and the fluid. 

To lowest order in the ratio a/R between bead radii 
and arm lengths (Stokeslet approximation [10]), the con- 
tribution to the trimer dynamics from the torque on the 
beads by the flow is disregarded. Faxen corrections at 
the bead scale and higher order images from interaction 
of the flow perturbation with the beads are disregarded 
as well. In this approximation, the equation of motion 
for the trimer can be written in the form 



X, = U(XO + U.(i) + F,/a, 



(2) 



where X^ is the coordinate of the i-th bead, Fj is the 
force on bead i by the rest of the trimer, a = dirfxa, with 
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H the fluid viscosity, is the Stokes drag, and Uj(i) is the 
flow perturbation generated in by the other beads. 

Contrary to the problem of the swimmer in a quiescent 
fluid, we see that the center of mass of the trimer X = 
(Xi+X2+X3)/3 plays a special role, which allows, in the 
language of [11], to fix the gauge for the problem. In fact, 
absence of external forces = 0, and linearity of the 

shear imply, from Eq. (2): X = U(X) + J2i U»(i)/3. In 
the absence of flow perturbation, the trimer would move 
as a point tracer located at its center of mass. Averaging 
over time the contribution by the flow perturbation gives 
the migration velocity of the trimer: 



A^2 



Ijmigr ^ (1/3)^(U.). 



(3) 



Similar conclusions can be drawn as regards rotation 
around an axis passing through the swimmer center of 
mass. The rotation frequency can be obtained from 
Eq. (2) imposing absence of external torques: X^ x 
F, = 0, that gives = /"^E^X, x [U(X,) + U,(i)], 
/ = Ei l-^iP- Neglecting the effect of the flow perturba- 
tion U, we flnd the constant rotation rate [12]: 



= a/2, 



(4) 



that coincides with the flow vorticity (the superscript in- 
dicates order in the deformation). This a bonus from the 
point of view of calculation that would not be available 
for less symmetric shapes (for instance, in the Stokeslet 
approximation, a dimer would tend to align with the 
flow). 

It is convenient to introduce a reference frame with 
origin in the center of mass, rotating solidly with the 

trimer, as illustrated in Fig. 1. Equation (4) allows to 
convert time averages into angular averages: 



{h) 



_ 1 



1 + -n^'\^) + . . . h{(i>), (5) 



with (f) the angle between the two reference frames (see 
Fig. 1). We shall use small case to identify vectors in the 
rotating frame. Equation (2) will become in this frame 
an evolution equation for the trimer internal degrees of 
freedom. 

The flow perturbation produced by the trimer can be 
expressed in function of the forces on the beads, and, 
thanks to linearity of low Reynolds number hydrodynam- 
ics, the resulting relation will be linear as well; in the 
rotating reference frame: 



(6) 



where T is the off-diagonal part of the Oseen tensor, re- 
lating, in the absence of an external flow, the Xj's and 
the fj's [10]. To lowest order in a/R, we can then write 
T(xj,Xj) ~ T(xjj), Xij = Xj - Xj, where, for i ^ j [10]: 



T(xij) = 



3a 
la 



(7) 




FIG. 2: Deformations of the trimer corresponding to Zi > 
for i = 1,2, 3. In the three cases Zj = for j ^ i. 



Substituting Eq. (6) into Eq. (3) and assuming exchange 
symmetry between the beads in the trimer: 



(8) 



where Ti = T(x2i) + T(x3i) and R, with Ru = R22 = 
coscj), R21 = —R12 = sin(^, is the rotation matrix back 
to the flxed frame. 

It is easy to see that a rigid trimer, due to the symme- 
try of its shape and of the external flow, would be unable 
to migrate. We must consider deformations. 

We can express the trimer deformation in terms of 
three independent parameters zi,Z2, Z3 as follows: 



/R 
/R 



{{V3/2){z2 + Z3),{zs-Z2)/2,0), 

{-{V3/2)z3,-zi-Z3/2,0), (9) 

(-(^3/2)^2,^1 +-22/2,0). 



As illustrated in Fig. 2, positive Zi corresponds to 
stretching of the arm opposite to bead i. From Eq. (8), to 
determine XJ™*^'", we must calculate fi at least to 0{z), 
z = zf)'^/'^ , but we can neglect the effect of the field 
perturbation. The rotating frame version of Eq. (2) reads 
at 0{(a/Rf): fi/a = x^ - u(xi). 

At 0(z°) there are no deformations and x*"' = 0. To 
next order, from Eq. (5), we can write x'^' = (a/2)x<^*' 
with prime indicating d/dct). From Eq. (4), the rotation 
velocity equals to 0(z*^) the vorticity of the field; hence 
u<°' is purely due to strain, u<'''(x) — Sx, with 

- - OL - - Oi 

Sii = -S22 = sin2(^, 512 = 821 = — cos2(j). (10) 

At 0{z), u receives a vorticity components from the 
correction to the trimer rotation velocity: u<^'(x) = 
— fi'^ X X. To 0{z) the trimer evolution equation will 
read therefore 



F(o) 



-aSx"', 



(11) 



f<^' = a[(a/2)xr''-Sxr'+f2'^'xxr']. (12) 

To obtain H'^', we apply x-"'x at both sides of Eq. (12) 
and sum over i, with the result 12*^'= Ei ^ Sx^^' / , 



3 



/<"> = J2t IxfT = -R^- Direct calculation from Eqs. (9- 

10) gives then 

Qfi' = |^\/3 (2:2-23) sin 2^+(2;2+2;3-22;i) cos 2^ ^. (13) 

Substituting Eqs. (5,7) and (9-13) into Eq. (8), we 
find the lowest order contribution V'^vr ^ (R(T^'"f{'' + 
-[-a)f^o)-j^(o)^ where indicates the uniform part of the 
average in Eq. (5). Imposing exchange symmetry among 
the beads: 



Zi = 'Y^[An COS iKpi + B„ sinn^i], 



(14) 



where = (t>, (62.3 = (f> ^ 2tt/3. We see that 
Eqs. (5,10,13) bring into Eq. (8) factors cos 2^ and 
sin 20 that combined with the coscp and sin^ in R, re- 
quire, for a non zero migration velocity, contributions 
oc cos 0, cos 3(?i, sin sin 3(?i in the deformation. Direct 
calculation, using T^" ~ (3{7/2 - [13zi + 29(z2 + Z3)]/8}, 
Ti2 ^ r^2i ^ /3(V3/8)(Z2 - 23) and /3{5/2 + [2zi - 

31(2:2 + -Z3)]/}8, ^ = 3a/(4a-i?), together with Eq. (14), 
gives in fact the result: 



migr 



■aa 



256 



73^1 + I3B3 



(15) 



The component Si corresponds to the pulsation regime 
depicted in Fig. 3. From here, the swimming strategy 
could be optimized minimizing at fixed JJ™^^^ the de- 
formation amplitude z or some more appropriate cost 
function [13]. 

Notice that, of all the contributions to migration, per- 
haps, only the one coming from the first term to right 
hand side of Eq. (12) qualifies as real swimming. It is 
in fact the only direct contribution from deformation to 
bead movement with respect to the fluid. All the others 
could be loosely described as a sort of sailing, in which 
the swimmer shape adapts to catch the external flow. 

The way the swimmer actually swims, is illustrated in 
Fig. 3. We see that propulsion does not come from direct 
drag by tJ, rather, from drag by the perturbation U gen- 
erated by interaction between the trimer and the external 
flow. Imposing exchange symmetry among the beads, al- 
lows, on the same line of Eq. (8) , to focus on the Stokeslet 
field of a single bead, say bead 1. Let us consider the 
contribution Bi sin 01 . We see that the Stokeslet field 
generated by bead 1 in response to the strain component 
of U, has a positive or negative component at beads 2, 3, 
depending on whether O<(/><7ror7r<0<27r. Migra- 
tion is produced by the deformation induced symmetry 
breaking between the two orientations. 

The swimmer of Eq. (14) behaves as a small engine, 
with the three arms stretching and contracting in se- 
quence. From linearity in z of J/""^'^, however, migration 
in the symmetric case of Eq. (14), implies migration also 
in the presence of a single deformable link. The scal- 
lop theorem does not apply in the present situation: the 
swimmer orientation and position are determined jointly 
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FIG. 3: The swimming strategy. When = 7r/2 (left), arm 
23 contracts, which turns out to maximize the drag along Xi 
by bead I's Stokeslet field [continuous line; see Eq. (7)]. The 
opposite occurs for — —tx/2 (right). Dashed lines identify 
the strain component of U. 



by the past history of deformations and of the gradients 
of U experienced: even correcting for the trivial displace- 
ment by U, reversing the deformations only, would not 
be sufficient to bring the trimer back to its original po- 
sition. Similar "violations" of the scallop theorem have 
been observed in [14], where the role of external flow was 
played by the perturbation generated in the fluid by other 
swimmers. 

This situation has a counterpart in the fact that swim- 
ming in an external flow occurs already at 0{z), and 
not at 0{z'^). The swimmer achieves this, exploiting 
the external flow to change orientation between reversed 
strokes; in a quiescent fluid, the orientation change would 
require other stroke components, and this is responsible 
for O(z^) swimming efficiency [11]. 

It is actually possible for the swimmer to extract the 
energy needed for migration directly from the external 
flow. The average power exerted by the swimmer on the 
fluid \sW = 3(xi • fi). Now, from xf = 0, the flrst 
contribution is W"'^' = 3(x^'> • f{°')(°), but the x^^' com- 
ponents contributing to swimming, entering Eq. (15), 
are odd in 0, while fj*"', from Eq. (11), is even. The 
same occurs, from Eqs. (5) and (13), with (.)*^'. Assum- 
ing that x<^' receives contributions only from ^1^3, the 
lowest order contribution to the power will be 



W-^^ = 3[(xi' 



(1) 



fa)^(o) 



f(0,^(0,]. 



(16) 



we see that choosing x*^' in appropriate way, we could 
set W^'^^ = 0, compensating to 0{z'^) the energy lost 
in swimming. From Eq. (11), the only terms in x'^' 
contributing to W^'^' are a cos 20, sin 20. Substituting 
Eqs. (9,11,14) into Eq. (16) we obtain: 



5B1B3 + 36s: 



3^ 



(17) 



where the inequality accounts for the internal dissipation 
of the device [at this point the condition that A2 = 0{z^) 
becomes irrelevant] . Looking at Fig. 3 and focusing again 
on bead 1, we see that this deformation corresponds to 
arm 23 being stretched when it is parallel to the flow 
direction, and contracted when it is perpendicular. En- 
ergy extraction from the flow comes from the fact that 
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FIG. 4: Swimming from braking. Suppose X23 enters the 
contracting quadrants in a stretched state (a). The brake 
is relased and the arm contracts, tiU it enters the stretch- 
ing quadrants fully contracted (6). The brake remains off 
till the arm returns, now fully stretched, into the contracting 
quadrants (c). At this point, the brake is acted on and |x23 
remains stretched until, passing through (d), it returns to (a). 



the |x23| grows in the stretching quadrants of the exter- 
nal strain, and decreases in the contracting ones (see Fig. 
3). 

What mechanism could allow a three-sphere swimmer, 
such as the one described, to passively undergo the de- 
formation sequence depicted in Eqs. (15) and (17)? A 
possibility is illustrated in Fig. 4: a system of orientation 
dependent brakes, controlling the stretching and contrac- 
tion of the trimer arms in response to the external strain 
(a more quantitative description of this mechanism is pro- 
vided in [15]). This mechanism contains all the required 



aspects of conversion of external flow to deformation and 
then to swim. We notice the asymmetry between the 
two configurations (6) and (d), that is necessary to mi- 
gration. At the same time, comparing orientations (a) 
and (c) with (b) and (d), we see that arm 23 is in the 
stretched state when aligned with U, meaning that en- 
ergy is extracted from the external flow. 

All the behaviors that have been described, derive 
from geometrical properties of the swimmer and could 
be expected to remain valid away from the small defor- 
mation regime considered. More importantly, it should 
be possible to construct a continuous counterpart of the 
trimer, for which the migration velocity predictions of 
Eq. (15) remain qualitatively valid, albeit substitution 
between the typical size of the moving parts in the two 
situations: a — > i?. Such a continuous swimmer should 
have the properties imposed by Eqs. (15) and (17): an 
overall shape in the laboratory frame that is stretched 
along U and asymmetric; say an egg with the tip point- 
ing at positive X2. In other words, a tank-treading cell 
[16, 17] with a fixed orientation asymmetric shape. Un- 
fortunately, this would require an orientation sensitive 
control of the device structure (compare with Fig. 4), 
and the practical realization of such a device at the mi- 
croscale would be very difficult. 

Nevertheless, from a conceptual point of view, some 
complex version of tank-treading appears to be the most 
effective strategy for a swimmer to exploit the presence 
of an external flow. A linear swimmer, say, could very 
well extract energy from U in a tumbling cycle, and use 
it to generate swimming strokes as if in a quiescent fluid. 
However, beyond having a lower stroke frequency than a 
comparable triangular device (it would spend most of the 
time aligned with the flow) , it would have only quadratic 
efficiency in the stroke amplitude. 
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